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By calculating orbitally resolved Pauli susceptibilities within maximally localized Wannier orbital 
basis transformed from first principles band structures, we find that magnetism in Fei+^Te still has 
its itinerant origin even without Fermi surface nesting, provided orbital modulation of particle-hole 
excitations are considered. This leads to strong magnetic instabilities at wave vector (0,7r)/(7r,0) in 
dxz/dyz orbitals that are responsible for the bicollinear antiferromagnetic state as extra electrons 
donated from excess Fe are considered. Magnetic exchange coupling between excess Fe and in-plane 
Fe further stabilizes the bicollinear antiferromagnetic order. Our results reveal that magnetism and 
superconductivity in iron chalcogenides may have different orbital origin, as Pauli susceptibilities 
of different orbitals evolve differently as a function of concentration of excess Fe and height of the 
chalcogen atom measured from the iron plane. 

PACS numbers: 74.70.Xa,71. 20.-b,75.25.-j, 74.25. Jb 



I. INTRODUCTION 

Unconventional superconductivity emerges from sup- 
pression of the magnetically ordered state in most iron- 
based super conductor J^IJ^. Thus, unveiling the origin of 
magnetism is a crucial step towards understanding the 
mechanism of superconductivity in these compoundJ^. 
In contrast to iron pnictides where magnetism can be ex- 
plained from both itinerant^ and localizecf^ limit, mag- 
netism in iron chalcogenides can only be accounted for 
by a localized scenario^ due to the fact that no Fermi 
surface nesting at the wave vector of (tt, 0, tt) is detected 
from angle resolved photoemission (ARPES) stud^EUll^ 
which conflicts with the observation of bicollinear anti- 
ferromagnetic (BAF) order from neutron diffraction anal- 
yg^ ji3 | i4j "Yhis remarkable difference between iron pnic- 
tides and iron chalcogenides is a serious challenge against 
establishing a unified theory for different families of iron- 
based super conduct or 

Han and Savrasov have attempted to recover the itiner- 
ant scenario for magnetism in Fei+^^Te by assuming that 
each excess Fe contributes eight extra electrons into the 
Fe plane, leading to a substantial shifting of Fermi level 
and a dramatic change of Fermi surface topolog}^^"^. 
A strong peak for condensation of particle-hole excita- 
tions emerges in the calculated Pauli susceptibility right 
at wave vector (tt, 0, tt) which accounts for the itinerant 
nature of BAF order. However, t his do ping effect was not 
observed in ARPES experiment ^1^^ and the oxidation 
state of interstitial Fe should be close to that of in-plane 
Fe as implied by neutron diffraction^^. Furthermore, our 
density functional theory (DFT) calculations based on 
nonmagnetic and BAF states both support the fact that 
all the irons share a similar valence state (see appendix 
[A[ ). Therefore the shift of Fermi level due to the extra 
electrons donated from excess Fe is one order of magni- 
tude smaller than that done by Han and Savrasov^*^ 



(see appendix B|. As a consequence, Fermi surface topol- 
ogy remains almost unchanged compared to the undoped 
case, leading to a good agreement with ARPES obser- 
vations. However, the calculated Pauli susceptibilities 
within a constant matrix element approximation as usu- 
ally done in the literatur d^ l ^^ l ^^ * only show a notable peak 
at (tt, tt, tt) in both undoped and doped cases, which again 
casts doubts on the itinerant picture of magnetism in 
Fei+a^Te (see appendices [C] and [D|) 

Here we will demonstrate that the ignored matrix ele- 
ments in calculated Pauli susceptibility are the key quan- 
tities for understanding the origin of itinerant magnetism 
in Fei+a^Te. The particle- hole excitation is strongly mod- 
ulated by orbitals, and its condensation wave vector is 
orbitally dependent. Within a reasonable amount of ex- 
tra electrons donated from excess Fe, prominent mag- 
netic instability at (0,7r)/(7r,0) in dxz/dyz orbital is ob- 
tained. The phase transition from commensurate to in- 
commensurate antiferromagnetic state observed experi- 
mentally can also be naturally explained as excess Fe 
further increases. Our results reveal that multiple insta- 
bilities coexist in different orbitals and evolve differently 
as a function of extra electrons and height of chalcogen 
atom measured from the iron plane, suggesting that mag- 
netism and superconductivity in iron chalcogenides may 
have different orbital origin. 



II. DETAILS OF OUR CALCULATIONS 

We use the experimental lattice structure of the para- 
magnetic phase^^ throughout the paper, except when we 
study the substitution effect of Te by Se, which can be ef- 
fectively viewed as a reduction of the height of the chalco- 
gen atom measured from the Fe plane. We employ the 
full potential linearized augmented plane wave method 
as implemented in WIEN2K^^ to calculate the electronic 
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FIG. 1: (Color online) Orbitally resolved Pauli susceptibility and spin susceptibility. Magnitudes of all the elements in Pauli 
susceptibility (a) and spin susceptibility (d) at g = (7r,0,7r). Dominant q-dependent Pauli susceptibility Xo^'^^(^5^ — 0) (b), 
Xo^'^^(^,C(; = 0) (c) at qz = tt. Dominant q-dependent spin susceptibility X^s^rpa(^^^ ~ ^) (^)' xt^RPAi^^^ ~ ^) (^) ~ ^• 

The two-dimensional contour maps are on the bottom. Here Fermi level is shifted up by 0.075eV, corresponding to the Fei+^Te 
compound at x = O.IG^^. Spin susceptibility is calculated at J = [//4 and U — 0.8 eV. The diagonal elements of Pauli 
susceptibility relevant to magnetism are marked by arrows in (a). From left to right, those are x^'^^^q^uj — 0) with t from 1 to 
5. 



properties of Fei+^^Te. An effective d-p model o f 16 o r- 
bitals within maximally localized Wannier basi d^^ * ^^ is 
constructed in order to calculate orbitally resolved Pauli 
susceptibility and spin susceptibility. We use a three- 
dimensional grid of 128x128x128 k and q points in the 
whole Brillouin zone with a temperature smearing of 0.01 
eV. Long range hoppings are involved so that both band 
structure and density of states from effective model are 
perfectly consistent with those from first principles in the 
energy window of [-6 eV, 3 eV] (see appendices [H| and [l[) 
The exact Pauli susceptibilit}^ at the Fermi level 
Xo^'*^(^7^ = 0) (defined in appendix [e| is a four-index 
tensor with each index varying from 1 to 5 which rep- 
resents five 3d orbit als of Fe (orbital 1: d^^^, 2: dxz^ 
3: dyz^ 4: (1^2 _y2 , 5: dz2^ where x, y, z refer to those 
for the original unit cell). In Fig. [l] (a), a total of 625 
elements of the Pauli susceptibility at the Fermi level 
are presented ai q = (tt, 0, tt) in a doped case where the 
Fermi level is moved up by 0.075 eV. It is found that off- 
diagonal elements of Pauli susceptibility are negligibly 
small, compared to the diagonal {p = s^q = t) elements 
which represent intraorbital and interorbital particle-hole 
excitations. Moreover, the spin susceptibility, which de- 
termines magnetism, is only related to the elements of 
Pauli susceptibility with the indices of s = t. Due to 
the above two facts, only intraorbital particle-hole exci- 



tations, i.e., Xo^'^^(<77^ — 0) with p = s = q = as 
marked by arrows in Fig.[l](a), have to be taken into ac- 
count in order to discuss possible magnetism, while other 
diagonal elements are irrelevant to the magnetism. We 
have also analyzed Pauli susceptibility at different wave 
vectors and in different doping cases including the un- 
doped case and found that the above conclusion is valid 
in all cases. 



III. RESULTS AND DISCUSSIONS 

In Figs.[l](b) and (c), two dominant elements of Pauli 
susceptibility Xq'^^{Q->^ = 0) with t = 3,4 are shown. 
The Fermi level is again shifted up by 0.075 eV in or- 
der to account for the excess electrons contributed from 
interstitial Fe. While a pronounced peak remains close 
to the wave vector q = (7r,7r,7r) in the d^'^-y^ orbital, 
a well-defined strong peak appears at g = (tt, 0, tt) in 
dxz orbitals, which is responsible for the BAF order. 

22*22 I I 

Xo ' (^7 ^ = 0) (see appendix [F|) shows a mirror sym- 
metry to y^'^^{q-,oj = 0) with respect to qx = Qy, which 
exhibits a strong peak ai q = (0,7r,7r). (7r,0,7r) and 
(0, TT, tt) instabilities are also found in the dxy orbital, 
which is considerably smaller than those in the dxz and 
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FIG. 2: (Color online) Evolution of particle- hole exci- 
tations in Pauli susceptibility at different shifted Fermi 
energies, (a) Xo'''\q.^ = 0). (b) xt^'^\q.^ = 0). 
AEf = 0, 0.045, 0.065, 0.075, 0.085, 0.105, 0.12, 0.15 eV, 
corresponding to the Fei+^Te compounds with x = 
0, 0.06, 0.09, 0.10, 0.12, 0.14, 0.16, 0.20^^. 

dyz orbitals. Similarly, the instability at (tt, tt, tt) in the 
orbital is much weaker than that in the dx2-y2 or- 
bital. These may indicate that dxy and d^^ orbitals play 
minor roles in the formation of magnetism or supercon- 
ductivity. 

Although it is interesting to get the fundamental 
(tt, 0, tt) and (0, tt, tt) instability by a reasonable rigid 
band shift served as doping extra electrons to the Fe 
plane from interstitial Fe, one may still doubt whether 
this instability vanishes as local Coulomb interactions 
and Hund's rule couplings turn on. In Figs, [l] (d), (e) 
and (f), we show the effect of interactions on the orbitally 
dependent instabilities. The spin susceptibility is a 5x5 
matrix and can be obtained in the form of Dyson-type 
equation within random phase approximation (RPA). (see 
appendix [E| Fig. [l] (d) presents all 25 elements of spin 
susceptibility at J/U = 0.25 with U = 0.8eV and at 
q = (7r,0,7r). It is found that off-diagonal elements are 
negligibly small, compared to the diagonal ones. As the 
J/U ratio decreases, while diagonal elements remain al- 
most unchanged, off-diagonal ones decrease and finally 
vanish at J/U = 0. The situation is the same for other 
wave vectors. Therefore, diagonal spin susceptibilities 




FIG. 3: (Color online) Evolution of particle- hole excitations 
in Pauli susceptibility as a function of Te-height measured 
from the Fe plane, (a) '^^(^,^ = 0). (b) Xo^'^^(^,^ = 0). 
Te height from hre = l.SA to 1.52A in the interval of 0.04A. 

play dominant roles in determining magnetism. Fig.[l](e) 
and (f) show two diagonal q-dependent spin susceptibili- 

33*33 44*44 

ties Xs rpa(^^^ ^ 0) rpa(^^^ ^ ^) Fermi 

level. We find that pronounced peaks are preserved at 
the same wave vectors in respective orbitals as what is 
detected in orbitally resolved Pauli susceptibility. We 
have also checked results at different values of U and 
found that the situation is the same. Therefore, we con- 
clude that for reasonable strength of interactions^^ the 
magnetic instability in spin susceptibility remains at the 
same wave vector as the condensation of particle-hole ex- 
citations in Pauli susceptibility. 

Next, we study the evolution of particle- hole excita- 
tions in Pauli susceptibility as a function of shifted Fermi 
energy. Fig. |2] shows momentum dependent Pauli sus- 
ceptibility of dxz and da.2_^2 orbitals along the path of 
(0,0, tt) (7r,0,7r) (7r,7r,7r) (0,0, tt) at different 
Fermi energy shifts, corresponding to different numbers 
of excess Fe. Without shift or as the shifts are small, like 
/\Ef = and 0.045, both xf '^^ and xt^'^^ show dom- 
inant peaks around (7r,7r,7r). While the peak in Xo^'^^ 
is slight enhanced at AEp = 0.045, compared to the 
case of AEf = 0, that in Xo^'^^ is strongly suppressed. 
At AEf = 0.065, 0.075, 0.085, a well-defined peak cen- 



4 



tered at (7r,0,7r) appears in Xo ' ^ indicating that the 
BAF state has its itinerant origin. The peak at (tt, tt, tt) 
in however, remains almost unchanged at these 

AEf. Further increasing A^;^ to 0.105,0.12,0.15, the 
peak around (tt, tt, tt) in Xo^'^^ rapidly decreases, while 
the peak previously right at (tt, 0, tt) in Xo ' nioves to 
an incommensurate wave vector, which is consistent with 
the experimental finding of transitions from BAF states 
to inc ommensur a t e phases as the number of excess Fe in- 
creases^^^'^^'^^'^. The larger the shift is, the farther away 
the peak is from (tt, 0, tt), which is also consistent with ex- 
periments. The incommensurate wave vector at (tt, e, tt) 
indicates that magnetic moments in each row along the 
a axis are antiparallel to each other, which modulates 
with propagating vector 27re/6 along the b axis. Here 
we pointed out that experimentally observed {6, 0, tt) in- 
commensurability may be related to the peak at incom- 
mensurate wave vector in Xo^'^^- We also check the tem- 
perature effect which can be effectively viewed as lift of 
the Fermi level. As expected, the peak moves away from 
(tt, 0,7r) as temperature increases, which also agrees with 
the experimental resultsf^. 

TABLE I: Comparison of energies among different antiferro- 
magnetic states. The coUinear antiferromagnetic (CAF), pla- 
quette antiferromagnetic (PAF), bicoUinear antiferromagnetic 
(BAF) and Neel ordered antiferromagnetic (NAF) states are 
taken into account within a supercell. Both local spin density 
approximation (LDA) and spin polarized generalized gradi- 
ent approximation (GGA) are used. FeTe denotes a supercell 
of FeieTeie. Fei.o625Te means FeirTeie while Zno.o625FeTe 
indicates ZnFeieTeie- (see appendix [gJ 

EcAF — EnAF EpAF — Enaf Ebaf — Enaf 

FeTe(LDA) -136.01 -111.84 -121.25 
FeTe(GGA) -138.22 -144.51 -162.83 

Fei .0625 Te 

(LDA) -126.41 -114.24 -135.81 

Fei.o625Te 

(GGA) -130.13 -151.00 -175.84 

Zno. 0625 FeTe 

(LDA) -128.79 -96.91 -107.19 



In spite of good agreements with various experiments, 
our itinerant picture of magnetism in Fei+^^Te still en- 
counters a severe problem. That is why (tt, 0, tt) in- 
stability wins the competition with (7r,7r,7r) instability 
which is obviously stronger than its rival. The reason 
is that excess Fe not only contributes extra electrons to 
the Fe plane which induces (tt, 0, tt) instability, but also 
provides a magnetic ion strongly coupled with in-plane 
Fe which stabilizes the BAF state, rather than the CAF 
one^^. In table 1, we show the importance of the in- 
terstitial magnetic ion. We construct a supercell with 
16 in-plane Fe. With such a supercell, various mag- 
netic orders like Neel antiferromagnetic (NAF), collinear 
antiferromagnetic (CAF), BAF, and plaquette antiferro- 
magnetic (PAF) order, can be studied on the same foot. 



(see appendix [G| for the cartoons) Without excess Fe, the 
magnetic ground state is strongly dependent on the func- 
tional. Local density approximation (LDA) yields a CAF 
state while generalized gradient approximation (GGA) 
favors a BAF state. This inconsistency casts doubt on 
the reliability of conclusions from previous DFT inves- 
tigations on FeTe where only GGA is used^*^. After 
putting one excess Fe into interstitial with the position 
according to neutron diffraction experiment J^^, we find 
that the BAF state becomes the ground state irrespective 
of which approximation one chooses, indicating robust- 
ness of the BAF state after involving excess Fe. The cru- 
cial role of existing magnetic ion in the interstitial can 
be verified by changing the interstitial Fe to Zn which 
contributes extra electrons with nonmagnetic ions. We 
find that the magnetic ground state becomes CAF within 
LDA. This is clear evidence that the existence of an in- 
terstitial magnetic ion is crucial for stabilizing the BAF 
state. On the other hand, (tt, 0, tt) instability also indi- 
cates a tendency towards the PAF state in principle, in 
addition to the BAF state. However, our results do not 
support the existence of the PAF state. 

Finally, we investigate the evolution of particle-hole ex- 
citations in Pauli susceptibility as a function of Te height 
measured from the Fe plane. We fix the number of ex- 
tra electrons to be 0.2e/Fe and tune the Te height from 
hre = I.SA to 1.52 A in the interval of 0.04A. From 
Fig. [3] (a), we find that lowering Te height rapidly sup- 
presses the (tt, 0, tt) instability in the dxz orbital. The 
peak position begins to move away from (tt, 0, tt) towards 
(J, 0,7r) at Htc = I.72A. The Xo^'^^ becomes featureless 
with further decreasing Te height, suggesting vanish of 
magnetic instability, which is consistent with the exper- 
iments where magnetic order disappears with substitu- 
tion of Te by Se^^^^. Note that since Se height is much 
lower than Te height, the instability to the BAF state 
is suppressed. Meanwhile, reduction of Te height also 
gradually suppresses the (tt, tt, tt) instability in the da,2_^2 
orbital as shown in Fig. Is] (b). However, it remains rel- 
atively large either at or close to (7r,7r,7r). Therefore, 
as the (7r,0,7r) instability of Xo^'^^ is completely sup- 
pressed, the remaining (tt, tt, tt) instability of Xo^'^^ ^^^Y 
be the source of superconductivity observed experimen- 
tall}'^, since its pairing is also associated with the real 
part of Pauli susceptibility through an integral over the 
Fermi surface. 



IV. CONCLUSIONS 

In summary, our results offer a possible explanation 
of why the Fermi surface is nested at (tt, tt) while mag- 
netic order at (tt, 0) in Fei+a^Te, as a hidden (tt, 0) in- 
stability coexists with the nesting of the Fermi surface 
at (tt, tt) within a reasonable range of extra electrons do- 
nated from excess Fe and height of Te ion measured from 
the Fe plane. Existence of interstitial magnetic Fe fur- 
ther selects the (tt, 0) instability as the itinerant origin of 
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the BAF state. Increasing excess Fe and decreasing the 
height of the chalcogen atom both suppress the (tt, 0) 
instabihty and may favor, respectively, the incommen- 
surate AF states and the superconducting states with 
pairing mediated by (7r,7r) magnetic fluctuations. Even 
though there exists a diversity of magnetic orders in iron 
pnictides and iron chalcogenides, their origins can be ex- 
plained from both itinerant and localized limit, suggest- 
ing that a unified model for describing iron-based super- 
conductors should involve both itinerant electrons and 
local spins. On the other hand, orbital differentiation has 
to be seriously taken into account, since different orbit als 
play different roles in magnetism and superconductivity 
as seen in Fei+^^Te where evolutions of Pauli suscepti- 
bilities of dxz/dyz and dx2_y2 are remarkably different 
when filling or structure is changed. Our findings may 
also point out a way of understanding magnetism and 
superconductivity in other iron-based superconductors, 
besides the existing localized scenario^^. 
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FIG. 4: (Color online)Possible valence of interstitial 

iron, (a) Partial density of state in nonmagnetic state. Solid 
line denotes the contributions from interstitial iron while dot- 
ted line the contributions from all the in-plane iron divided 
by a factor of 2. (b) Partial density of state in bicoUinear 
antiferromagnetic state. Now the contributions from all the 
in-plane iron are divided by a factor of 4. The supercell we 
used is FeirTeie. 



Appendix B: Shift of Fermi level VS number of 
excess iron 
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Appendix A: Possible valence of interstitial iron 

In order to estimate possible valence of interstitial 
Fe, we analyze partial density of states obtained from 
both nonmagnetic state and bicollinear antiferromagnetic 
state of a supercell of FeirSeie with sixteen in-plane Fe 
and one interstitial Fe. The supercell is constructed ac- 
cording to experimental structure of Fei.oesTe at 80K!^. 
The results from first principles calculations^^ are pre- 
sented in Fig. [4j We find that 3d orbitals of interstitial 
Fe are also prominently occupied, indicating that inter- 
stitial Fe should not be in a valence state of +8 where 
3d orbitals of interstitial Fe should be empty. Instead, 
we find that the occupation number on interstitial Fe is 
almost the same as that on in-plane Fe by comparing 
the integrated density of state under Fermi level for all 
the irons. Therefore, possible valence of interstitial iron 
should be the same as that of in-plane Fe, i.e. ~ +2. 
Moreover, the density of state of interstitial Fe in non- 
magnetic state (Fig. |4](a)) shows a strong peak at Fermi 
level, indicating a strong tendency towards a magneti- 
cally ordered state. 



Based on the fact that valence of interstitial Fe is +2, 
extra electrons doped into Fe plane should be counted 
as 2x in Fei+a^Te compounds, rather than sd^^''^. By 
rigid band shifts, in Fig. [5j we show shift of Fermi level 
as a function of concentration of excess irons. The lat- 
tice constants are taken from experiments^^. The height 
of Te ion measured from Fe plane is varied from 1.52 A 
to 1.80 A in a interval of 0.04 A, since it is well known 
that the electronic properties of Fe-based superconduc- 
tors are strongly dominated by the height of anion. By 
applying first principles calculations, we find from Fig. [5] 
that the shift of Fermi level should be one order of mag- 
nitude smaller than what has been done by Han and 
Savraso^PHl^I. 



Appendix C: Fermi surface 

In Fig. [6l we compare the Fermi surfaces in the ab- 
sence and presence of Fermi energy shift. We find that 
with a reasonable Fermi energy shift, such as 0.075eV, 
while the radius of Fermi cylinders at the F and M point 
are different, the overall shapes of Fermi surfaces are still 
the same, indicating that Fermi surface topology is still 
consistent with that observed by angular resolved pho- 
toemission spectroscopy even though there is a shift of 
Fermi energy due to the extra electrons contributed from 
interstitial excess Fe. 



Appendix D: Pauli susceptibility within Constant 
matrix element approximation 

In Fig. [7| we show the results of Pauli susceptibilities 
within constant matrix element approximation. Without 
a Fermi level shift, a strong peak at q = (7r,7r,7r) is ob- 
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FIG. 5: (Color online)Shift of Fermi level as a function of 
excess iron concentration. Experimental lattice constants 
are used while height of Te ion measured from Fe plane is 
varied from 1.52 A to 1.80 A in a interval of 0.04 A. The 
Shift of Fermi level of Fei+a^Te compound at x = 0.068 calcu- 
lated from density functional theory is one order of magnitude 
smaller than what has been done by Han and Savraso\ff^^. 




FIG. 7: (Color online)Pauli susceptibilities. Here we adopt 
constant matrix element approximation, (a) Without shift of 
Fermi level, (b) Fermi level is shifted up by 0.075eV. Exper- 
imental lattice structure at high temperature is used^. 



tight-binding Hamiltonian. Here f{E) is the Fermi dis- 
tribution function. The spin susceptibilit}^^^ within ran- 
dom phase approximation (RPA) is obtained in the form 
of Dyson type equations as 



yCs,RPA 



Xs,RPA^s Xo l^^J 



where repeated indices are summed over. 




(JJlJ. 



FIG. 6: (Color online)Fermi surfaces. (a) Without shift of 
Fermi level, (b) Fermi level is shifted up by 0.075eV. Besides 
the difference in the radius of Fermi cylinders at the F and 
M points, the overall shapes of Fermi surfaces in these two 
cases are quite similar. Experimental lattice structure at high 
temperature is used^. 



served. With a Fermi level shift of 0.075eV, the peak at 
q = (7r,7r,7r) becomes a bit broader. Howeverno promi- 
nent peak can be observed dit q = (tt, 0, tt) or q = (0, tt, tt). 



Appendix E: exact expression of Pauli susceptibility 
and spin susceptibility within RPA 

The Pauli susceptibility^ is defined as 

^ ^ a^(/c)a^*(/c)<(/c + g)at*(/c + q) 



AO 



E,{k + q)-E^{k) + i()+ 
x[f{E,{k + q))-f{E^{k))] (El) 



N ^ w ■ 



where matrix elements a^(/c) = {s\fik) connect the orbital 
and the band spaces and are the components of the eigen- 
vectors resulting from the diagonalization of the effective 



Appendix F: Orbital resolved Pauli susceptibility 
and spin susceptibility 

In Figs.[8](a)-(c), three elements of Pauli susceptibility 
Xo^'^^(<7, cc; = 0) with t = 1, 2, 5 which is relevant to mag- 
netism are shown. In Figs, [s] (d)-(f), corresponding spin 
susceptibilities are present. 



Appendix G: Cartoons for different magnetic states 

In Fig. [9) we show the cartoons for different magnetic 
states we considered in this investigation. Please note, 
the (tt, 0) instability favors two kinds of magnetically or- 
dered states. One is the bicollinear antiferromagnetic 
(BAF) state (Fig. |9] (c)), the other is the plaquette an- 
tiferromagnetic (PAF) state (Fig. |9] (d)). From our first 
principles calculations, the energy of the BAF state is 
always lower than that of the PAF state in all the cases. 



Appendix H: comparison of Band structure and 
DOS from LDA and GGA 



Fig. [10] shows the comparisons of band structure, as 
well as density of state, obtained from local density ap- 
proximation (LDA) and generalized gradient approxima- 
tion (GGA). The results from these two approximation 
are perfectly consistent with each other, indicating that 
our results are independent of which approximation we 
choose. 
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FIG. 8: (Color online) Orbitally resolved Pauli sus- 
ceptibility and spin susceptibility, q-dependent Pauli 
susceptibilities xj^'^\<?,^ = 0) (a), xf'^^(<?,^ = 0) (b) 
and Xo^'^^{q^^ — 0) (c)- q-dependent spin susceptibilities 
xI]rpa(Q^^^ 0) (d), x1^rpa(Q^^^ 0) (e) and x^^ip^(^,cJ = 
0) (f). Here Fermi level is shifted up by 0.075eV. Spin sus- 
ceptibility is calculated at J — U /A and [/ = 0.8 eV. 
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FIG. 10: (Color online)Band structure and DOS from 

LDA and GGA. (a) Comparison of band structure from lo- 
cal density approximation (LDA) and generalized gradient 
approximation (GGA). (b) Comparison of DOS from LDA 
and GGA. 
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FIG. 11: (Color online)Band structure and DOS from 
DFT and effective model, (a) Comparison of band struc- 
ture calculated from density functional theory and effective 
tight banding model, (b) Comparison of DOS calculated from 
density functional theory and effective tight banding model. 
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FIG. 9: (Color online) Cartoons for different magnetic 
states. (a) Neel ordered antiferromagnetic (NAF) state, (b) 
CoUinear antiferromagnetic (CAF) state, (a) Bicollinear anti- 
ferromagnetic (BAF) state, (b) Plaquette antiferromagnetic 
(PAF) state. 



Appendix I: comparison of Band structure and DOS 
from DFT and effective model 



Fig. [TT] presents the comparisons of band structure, 
as well as density of state, obtained from first princi- 
ples calculations and tight-binding model calculations. 
The effective tight-binding model is derived from first 
principles band structure through a transformation from 
Bloch space to maximally localized Wannier orbital ba- 
sis. The energy window we used is from -6 eV to 3 eV 
where the weight of each eigenstate is mostly contributed 
from Fe 3d orbit als and Te 4p orbit als. In order to get 
a perfect consistency in band structure calculated from 
first principles and the effective tight-binding model, we 
keep the long-range hoppings from the atoms in the unit 
cell at [0,0,0] to those at [9a, 9b, 5c]. The results cal- 
culated from density functional theory and from the ef- 
fective tight-binding model are perfectly consistent with 
each other, indicating that the transformation we used 
will not impose any systematical error on our results. 
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